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In this work, we study the collapse dynamics of an inhomogeneous spherically symmetric star 
made of dark matter (DM) and dark energy (DE). The dark matter is taken in the form of a dust 
cloud while anisotropic fluid is chosen as the candidate for dark energy. It is investigated how dark 
energy modifies the collapsing process and is examined whether dark energy has any effect on the 
Cosmic Censorship Conjecture. The collapsing star is assumed to be of finite radius and the space 
time is divided into three distinct regions E and , where E represents the boundary of the star 
and V~ {V^) denotes the interior (exterior) of the star. The junction conditions for matching 
\^ • over S are specified. Role of Dark energy in the formation of apparent horizon is studied and central 

' singularity is analyzed. 

o : 

CN ■ PACS numbers: 04.20 CV, 04.70 BW, 04.20 DW. 

D ■ 

C/0 ■ I- INTRODUCTION 

<N ■ 

^-H ' The collapse dynamics in the framework of Einstein gravity was started long ago by Oppenheimer and 
Snyder [1]. An extensive study of gravitational collapse for matter cloud in the form of pressure- less dust 
• has been done in recent years [2-6] (for recent reviews see [7]). These studies show the validity of Cosmic 
' Censorship Conjecture in six and higher dimensions under reasonable physical conditions, the role of initial 
data on the end state of collapse and some dynamical symmetries [8] on the initial data set. Also the collapse 
'[ dynamics has considerable astrophysical significance. 

' So far, these studies are mostly concentrated to dust cloud and matter with pressure components [9-13]. 

Also there are few works [14-15] where in addition to dust matter, one has cosmological constant with it. 
' Recently, Madhav, Goswami and Joshi [16] have studied in details the effect of a cosmological constant in the 
Q I background of asymptotically anti de-Sitter (or de-Sitter) space -time. They have treated the cosmological 
O^" constant as a dark energy component, motivated by the recent astronomical observations [17] of high redshift 
I type la Supernova. Usually, the dark energy behaves as the source of repulsive gravity and is considered to be 
bJO important in the present era to justify the undergoing accelerated expansion of the universe. 



X 



Thus, it is interesting to study gravitational collapse of a matter cloud containing dark energy. In the present 
work, we examine the collapsing process of an inhomogeneous spherically symmetric star whose matter inside 
^ , contains a combination of dark matter (dust) and dark energy. We have studied in two different sections 
- - ' the collapsing scenario for TBL model and general spherically symmetric model. The formation of apparent 
horizon has been investigated in details for various cases and it is examined whether bouncing solutions are 
possible or not. The paper ends with a short discussion and concluding remarks. 



II. BASIC EQUATIONS 

The metric ansatz of a general spherically symmetric metric can be written as 

ds^ ^ -e^'^dt^ + e^'^dr^ + R^dil^ (1) 

where a, 7 and R are functions of t and r and dil^ — d9^ A- siv? 9d<f>^ is the usual line element on a unit 
2-sphere. The form of the energy momentum tensor is taken as 
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A* (total) " A* (DM) ^ M (DE) ^ ' 

Here the coordinates {t, r, 6, (f) are taken to be comoving and the components of the above energy momentum 
tensor in this comoving frame are 

(DM) =^^«5(-PoM, 0,0,0) (3) 

and 

'^!l(DE)='^^"'9i-P,Pr,Pt,Pt) (4) 

The quantities p^^ , p, pr and pt are respectively the dark matter density, dark energy matter density, radial 
and tangential pressure for dark energy and are functions of r, t. In order to satisfy the weak energy condition 
for the whole matter we have 

for any time- like vector V^. In explicit form this means 

Pt>0, p^+Pr>0, pT+Pt>0 (6) 

where p^ = + p, is the total matter density. Now the explicit form of the Einstein's field equations for 
the metric (1) with matter field given by (2) are (choosing SttG = c = 1) 



(8) 



-e-'Wa + d^-dT+f + <i|-7|) +e-»°(7" + 7"-aV + ^+7'^-c.'f 1 =» P) 



and 



R' - dR' -j'R = (10) 

where over dot(.) and dash(') stand for partial derivatives with respect to t and r respectively. We now 
introduce the mass function m{t,r), defined by 

m{t, r) = R{1 + e-^^te - e'^^R'^) (11) 

which can be interpreted as the total mass inside the comoving radius r at the instant t. This definition was 
first introduced by Cahill and McVittie [18] and since then it has been widely used in different context. Also, 
for asymptotically flat space times, it corresponds to the correct Bondi mass at infinity. If the boundary of 
the collapsing cloud is leveled by the comoving coordinate rs, then the total mass of the collapsing star at the 
instant t is given by 
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TOs(i) =m(i,rs) (12) 
For regularity on the initial hypersurface t — ti, we must choose ■m{ti,0) = 0. 

Using this mass term, the field equations (7) and (8) can be written in compact form as 

m' rh 

As the matter field is a mixture of dark matter (inhomogeneous dust) and dark energy (anisotropic fiuid), so 
from energy momentum conservation relation T^" = 0, we obtain 

f r\ R 
p^+p^ ia + 2-\ +dpr + 2pt-=0 (14) 

and 

R' 

Pr + {Pt + PrW - ^{Pt -Pr)^=0 (15) 

If Q{t, r) denotes the interaction between the dust cloud and dark energy, then a self consistent conservation 
equation for each of them can be separately written as (see eqn(14)) 



(16) 





and 



---Q{t,r) (17) 

Thus we have six independent differential equations namely equations (10), (11), (13), (16) and (17) 
containing eight unknowns (geometrical or physical) namely, p^^j^j , p,pr,pi, Q{t,r), a,"/ and R, giving us 
freedom to choose two free functions. In particular, for a given initial data satisfying the weak energy condition, 
the choice of these free functions will determine the matter distribution and metric of the space time and hence 
a particular dynamical evolution of the initial data. 

On the initial hypersurface: t — ti wc have to specify eight functions of r namely 
PDM„i'r)^Po{'r),Pr„{r),Pt„ir),jo{r),ao{r),Ro{r) and Qo{r). By proper scaling of the radial coordinate r, 
we can choose Ro{r) = r, without any loss of generality. It is to be noted that the remaining seven (initial 
data) functions arc not independent but are related by the field equations. Further, to preserve regularity and 
smoothness of the initial data, we have to assume that the initial pressures at the regular centre r = should 
satisfy 

i) gradient of the initial pressures vanish at the centre i.e, ^^(,(0) = ^(^(0) = 

ii) the initial pressure should be isotropic at the centre i.e, Pro(O) = Pto(O)- 

From the field equations (13), we see that the total matter density becomes infinity either when i? = or 
i?' = or both. The case for i?' = is termed as shell crossing singularity which is not of much physical 
interest as it is gravitationally weak in nature and is removable singularity. The other one (which we shall 
discuss here) is termed as shell focusing singularity where the physical radius of all matter shells goes to zero 
{R = 0). So we shall assume i?' > at all instant to avoid any crossing of the shells. If ts{r) denotes the 
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instant when a shell at coordinate radius r collapse to r = 0, then we have R{ts{r),r) = 0. ( Note that the 

singularity time tg is a function of r due to the inhomogcncity of the space-time ). In the subsequent sections 
we shall consider the collapse dynamics for two separate assumptions namely 

i) a geometrical assumption: 7 = (LTB model) and 

ii) a physical assumption: Pr = 0. 



III. COLLAPSING PROCESS IN LTB MODEL WITH DM AND DE 

The metric of the space time describing the collapsing matter in the standard Lemaitre-Tolman-Bondi (LTB) 
model is given by 

ds^ = + ^'l, , dr^ + R^dn^ (18) 
1 + j{r) 

The simplified field equations are 

7?2 RR^_J_ f__ 

R"^^ RR'~ R^~ RR'~^^ ^ ' 



i?2 R2 R 

and 



f R^ „R /„„^ 



'2RR^ R R' R R 

with energy conservation equations 

' .R k' 



f R H! R R^ 

7=Pt (21) 



Pbm+ 2- + — U^„=Q(i,r) (22) 



I ^R R' \ R' ^ R ^ / X /^^x 



and 

p'r = '^{Pt-Pr)^ (24) 

From equation (24), it is to be noted that 

a) anisotropy of pressure implies the inhomogeneity of the pressure and vice- versa 
b) dark energy model with vanishing radial pressure and non-zero tangential pressure [pr = 0,Pt 0) is not 
possible. 

From the recent study of the role of pressure in quasi-spherical gravitational collapse, we choose the radial 
pressure as [12] 
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Pr(t,r)=-^ (25) 



Then from the conservation equation (24), the tangential pressure has the expression 



^ (+ ^\ sir) /, n\ g'{r) 

(Note that we have taken — ve sign for pr as pressure is — ve for dark energy). The arbitrary function g{r) 
should be Qir"^"^) near r = for the regularity and smoothness of the initial pressure. Solving the field 
equations (13) we have the expression for mass function and total matter density as 

m{t,r)=mo{r) + j^^^y^, {n^3) (27) 

and 

TOn(r) Q'(r) q(r) , ^ 

^^(^'^) = W+ (3^1)^ + 1^ (28) 

where mo(r) is an arbitrary function of r, restricted by the energy conditions. The evolution equation for the 
area radius R can be obtained from the field equation (20) (or from the definition of the mass function (11)) as 

Now using the scaling independence of the coordinate r [16] let us write 

R{t, r) = rv{t, r) (30) 

so that 

v{ti,r) = l, v%,r) = Q, v <Q (31) 

Then from (29) the time evolution equation for v is 



/v=^ (32) 



(3 - n)r"-i 

This is the evolution equation of a particular shell and V{r, v) may be termed as an effective potential. 

One of the main questions that we shall address here is how the collapse dynamics gets modified due to 
the presence of dark energy. In case of pure dust collapse, once the collapse initiates, starting from an initial 
regular hypersurface, the cloud necessarily collapses to a singularity due to gravitational attraction without 
any reversal or bounce. On the other hand, the negative pressure of the dark energy produces repulsive 
gravitational force and slows down the collapsing process. So it is expected that, due to the presence of dark 
energy, the collapsing process may not be smooth. Thus, it is interesting to study in details the qualitative 
nature of the effective potential V{r,v) for which the allowed regions of motion correspond to V{r,v) > 0. In 
fact initially, when the collapse starts, we have w < and we will have a rebounce if we get v = Q before the 
shell has become singular. In other words, study of various evolution for a particular shell is equivalent to 
study the nature of the roots of the equation V{r, u) = for fixed r. 
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To get a clear idea, let us consider a smooth initial data for which initial density, radial pressure and curvature 
are smooth near r = and the corresponding functions will have power series expansion as follows: 



Pto {ti, r) = Poo + po2r^ + Po4r* + 



Pro{ti,r) = -poJ^'^^'r'^ 
f{r) = hr^ + hr^ + ■ 



Poi^'^^'r^ -po(i^''V 



(33) 



As we are concentrating on the evolution of the shells near r = so we may neglect here higher order terms 
in the above expansions. Thus near the centre (r << rs) equation (32) becomes 



/v (34) 

The numerator of the right hand side of the above equation (i.e the effective potential V(r,v)) is a 
polynomial in v and in general may or may not have positive real roots which correspond to physical cases. 
As V{r,{)) = ^ > so any region between R = rv = and the first positive zero of V{r,v) always becomes 
singular during collapse while if there are two consecutive positive real roots then the region between them is 
forbidden as i;^ < there. Hence a particular shell will bounce if it lies in the initial epoch {v = 1) to the right 
of the second positive root. We shall present various possibilities for n = 1 and derive the necessary conditions 
below: 

a) /o = 0: (Marginally Bound Case) 

Here as both poo and po2^^'^ are positive so both roots of the quadratic equation are complex and hence v'^ is 
always positive and no bounce occurs. Hence a singularity always forms from the initial collapse. 

b) /o > 0: 

Both roots will be real and negative if fo^ > ^poo Po2^^^'r''^, otherwise they are complex in nature. 

c) /o < 0: 

For the restriction fo^ > |poo Po2^^^r'^, both roots are real and positive, otherwise they are complex. So 
bounce will be possible in this case. 

Horizons: 

For observers at infinity, the event horizon plays an important role in characterizing the nature of the singu- 
larity. As the calculations are complex in nature, so instead of event horizon, one considers a trapped surface 
which is a compact space-time 2-surface with normals on both sides as future pointing converging null geodesic 
families. In fact, if (r =constant, t =constant) the 2-surface Sr.t is a trapped surface then it and its entire 
future development lie behind the event horizon provided the density falls off fast enough at infinity. As appar- 
ent horizons (AH) are the boundaries of trapped regions (surfaces), so in the present case it can be written as [19] 

AH : g'^-'R^R^ = 

Using the metric form (1) one gets 

which with the help of the mass function simplifies to 



fov + — + — -V' 

3 (S — n) 



^ m{t,r) _^ 



(35) 



In the present case (i.e, for LTB model ) the apparent horizon is characterized by 



(36) 



(37) 



We shall now discuss the formation of apparent horizon for different values of n: 
n = l: 

In this case, equation (36) is a quadratic equation in R, having real positive roots if 
(i) < mog < \ and (ii) g{r) > 
and these horizons are given by 

i?e(r) = {1- Vl-2mo5}/5 
Rb{r) = {l + Vl-2mo5} /g 
These are usually termed as the cosmological and the black hole horizons [19] . 

n = 2: 

Equation (36) reduces to a linear equation in R, so there is only one horizon namely 
Rbc{r) = (assuming g{r) < 1). 

n = 4: 

Again equation (36) is a quadratic equation in R, having two horizons if |mo| > 2y^. 
n = 5: 

Here equation (36) simplifies to a cubic equation in R as 

2R^ - 2mo(r)i?2 + g(r) = (38) 
In the following we shall discuss the nature of the roots of the cubic equation for various possibilities: 

i)Forfl(r)<i^(l + ^), 

there are two positive roots of (38), which correspond to two apparent horizons namely the cosmological 
horizon and the black hole horizon given by 



_ mo(r) 2rno'' 
(27)3/ 



R,{r) = HMIi + 



■Cos 



2mo^ V 2 27 



}]J 



(39) 



ii) When5(r) = ^(l + ^), 



then there is only one positive root (corresponds to a single apparent horizon) of the above cubic equation 
given by 



^fec(r) = ^ + moV(27)i 



(40) 



iii) Ifg(.)>l^(l + ^), 

then there are no possible roots and hence there are no apparent horizons. 

Further, if tah{i~) be the time of formation of apparent horizon then from the evolution equation (34) we have 

- tsir) = - r ^ r (41) 



where we have considered shells close to r = 0. Here ts{r) is the time of formation of singularity of a shell at 
comoving radius r and Vah = Rb (or Rc)/r. 

It is to be noted that the presence of the dark energy in the form of anisotropic fluid characterizes the 
formation of apparent horizon and also influences the time difference between the formation of apparent horizon 
and singularity formation. 

The above time difference will characterize the final end state of gravitational collapse and also the nature 
of the resulting singularity. If the formation of the horizon precedes the formation of the central singularity, 
then the singularity will necessarily be covered, i.e, black hole will be formed. On the other hand, if the time 
difference is reversed then end state of collapse leads to a naked singularity. As this characterization of the 
singularity is local, so for global visibility, we should examine whether it is possible to have future directed null 
geodesies that terminate at the singularity in the past. Since analysis of null geodesic has been done widely so 
we are not presenting it here. 



IV. GENERAL SPHERICALLY SYMMETRIC COLLAPSING MODEL WITH = 0,pt ^ 

FOR DARK ENERGY 

Over the last few years spherically symmetric collapse with anisotropic pressure in the form of vanishing 
radial pressure have been studied in details. The main objective of the present study is to examine whether 
bouncing situation is possible due to the presence of dark energy and to investigate singularity formation 
conditions. Also we try to understand how dark energy affects the junction conditions at the boundary of the 
collapsing star to an exterior region, the formation of trapped surfaces (apparent horizon) and the nature of 
the central shell focusing singularity. 

For general spherically symmetric model (given by equation (1)), if we assume the radial pressure to be zero 
then from the field equations (13) we have 



..(M) = ^ (42) 



and from the conservation equation (15) using (42) we get 



Vt = (43) 

Here m(r) is an arbitrary function of r, satisfying the energy conditions. For the remaining freedom to choose 
one function, we take the metric coefficient 7(r, t) in the specific form 



7(r,i)=7o(t)+/3(i?) 
(motivation for choosing such a form is given in ref. [16]). 
Using this choice of 7 in equation (10) we have the other metric coefficient as 



(44) 
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^ R'it,r) 



with b{r) an arbitrary (7^ function of r. 



These choices of the metric coefficients when substituted in the definition of the mass function (i.e, equation 
(11)), we obtain the evolution equation for the area radius R as 



R' 



p2/3(fi) 

R 



m(r) + {62(r)e2'3(«) - l} 



R 



(46) 



where a multiplicative function of time has been eliminated by a suitable scaling of the time coordinate. As 
in the previous section, introducing R{t, r) = rv{t, r) the above equation transforms to 



V = 



m(r) {b^(r)e^^(''^) - l} ' 



(47) 



To study this evolution equation near the singularity (r = 0), one should consider smooth initial data (i.e, 
initial density and pressure) with power series expansion (near r = 0) 



Pi 
Pti 



{ti, r) = po+ P2r^ + Pir'^ + ■ 



■■ Pt{U,r) = -Pt2r^ - PtAT^ - Pter^ 



and 



b^r) = 1 + bo2{ry + boiiry + . 



(48) 



where the first term in the series expansion for 6^(r) is chosen to be unity as the metric looks Minkowskian 
near the centre r = 0. Thus close to r = 0, neglecting higher order terms in the above expression (47) for i), 
one gets 



•2 A 2pt2 2 2\ f Po , , '2pt2bir) 3\ 

V = [1 r^v^ Hr + bo2V u /v (49) 

V Po / V 3 Po J 

During the collapsing process the scale factor v decreases from unity (at the initial epoch) to zero (at the 
time of singularity), so the first factor in the numerator of the right hand side of (49) is initially positive (as 
it represents e^^). Hence this factor is positive definite during the collapse dynamics. Thus, whether a shell 
bounces or not is completely determined by the second factor in the eff'ective potential. 

The second factor, which is a cubic equation in v, has in general three distinct roots. As positive real roots 
are only physically admissible, so for different situations wc; cixaniine below whether positive roots are possible 
or not. As the term in quadratic power of v (i.e, v"^) is absent, so whenever all the three roots are real, there 
should be at least one positive and one negative root. 
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TABLE (Analysis of the roots of the equation) 



Condition 


Nature of the roots 


Physical Consequences 


(i) &02(r) > 


Exactly one possible root 
(say a{r)) and the two 
other roots are either negative 
or complex conjugate 


The physically allowed range of v is 
[0,a]. If a > 1, 
then there will always be a 
singularity while a < 1 implies 
an unphysical situation initially as 

v'^ < 0, i.c, all shells in 
the allowed dynamical space [0, a] 
becomes singular starting from initial collapse 


(ii) 602 (r) < 
(a) po > I """^1 

\ J 1^^ 3 Pt2b{r) 


No positive real root 


Singularity is always the final 
outcome of the collapse 


CM n "2 bo2^(r) 
(bj PO < 3 p^^b{r) 


Two positive roots vi (r) 
and V2{r) {vi < V2) where 
t;i = A Cosip, 
V2 = X Cosii; + 4f ) 

and 3^ = Cos-^ [- 


Allowed space for dynamical evolution 

is when V lies in the range [0, i-'i] 
and [v2,oo). The value of v in 
{vi,V2) is not allowed as 
v"^ < 0. Initially, shells for 
V £ [0,vi] always become singular 
but shells for v G 00) will undergo 
a bounce and subsequent expansion starts 
from initial collapse 


C„\ „ _ 2 bo2''*(i-) 
(Cj PO - 3 pt26(r) 


Two equal positive roots 


There is no forbidden region 
and bounce will occur if A < 2 



In case (b)(of (ii)), there are two positive real roots and bounce will occur when the area radius approaches 
= rv2. As in this case both singularity and bounce may occur so the explicit conditions for a particular 
shell to become singular or undergo a bounce can be stated as 

i) For singularity: A Costjj > 1 

ii) For bounce: A Cos(V' + ^) < 1 

So far, we have analyzed the evolution of shells close to the centre. But it is difficult to study (analytically) 
shells far from the centre and simple expression for the singularity or bounce is not possible, only one has to 
use numerical methods. However, the above analysis can be applicable to the entire cloud within the boundary 
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r = r-E, provided either the initial area radius rs is itself small compared to the initial data coefficients, i.e, 
Pi^s*j Pti^T.^ « 1 or the initial data are such that higher order coefficients in the power series expansion are 
identically zero i.e, pi = = P(i,V i > A and P2?'E^, Parn^ « 1- Moreover, if the whole collapsing cloud 
undergoes bounce starting from the initial collapse (i.e, case(b) of (ii) holds for all r in < r < rs) then the 
suflacient condition for avoiding shell crossing is ^^}^(^)^ > 1)V r G [0, Te) and (5 > 0, a small increment. 



V. SPACE-TIME MATCHING 

We shall now show the space-time matching of the collapsing star with the exterior vacuum. First of all, we 
shall find the matching conditions for a general exterior space-time and then for Schwarzschild space-time as a 
particular case. For smooth matching we shall use the Israel-Darmois [20] junction conditions on the boundary 

E- 

Let the interior of the star is denoted by V~ and is the exterior vacuum space-time region with the 
surface of the star as the common boundary. So the interior space-time {V~) is described by the metric given 
in equation (1) while that of the exterior vacuum {V~^) is taken as 

V+ : ds+^ = -A^{T, A) dT^ + B'^{T, A) [dh? + A^dfi^] (50) 

The equation of the bounding surface Y^, considering it as an embedding in the interior or exterior space-times 
are given by 

E- : r = rE; E+ : A = As(r) (51) 
with rs, a constant. Then the metric in ^ relative to the coordinates of V~ and are respectively 

S" : dss- ^ = e-^^df + R^{t, rs) dn^ (52) 

and 

E+ : = -A^{T, As) dT^ + B'^{T, As) [dAs^ + As^ dU^] (53) 

According to Israel-Darmois, the matching of the interior and exterior demands the continuity of the first 
and second fundamental forms on the bounding hyper surface. (Note that, if there are matters on the boundary 
then the second fundamental form has a jump discontinuity over the boundary): 

The continuity of the first fundamental form gives 

2 

dT = dr (54) 
For matching of the second fundamental form we need extrinsic curvature which is given by 

where C are the intrinsic coordinates on x'^ are the coordinates of V^, the ( , ) denotes partial diflFeren- 
tiation and the expression for the unit normal ricr to the hypersurface (/ =constant) is given by 



i) BAs = R{t,r^) = i?E (say) ii) e^dt = 



rfAs 
dT 



In the present matching model the explicit form of the unit normal is given by 



(56) 
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2,- = (0,e"^''--\0,0) and n<,+ = (^-^, 1, 0, 0^ 



AB 



Then the non- vanishing components of the extrinsic curvature in y~ and have the expressions 



(57) 



(58) 



V+ 



Ktt^ = 



AB 



B dB 
A^df 



dT 



2 dA 
AdX 



1 dB 
B'dA 



V dT 



and 



_^ 1 dA 2 dB 



AdT B dT J \ dT 



dT2 



(59) 



AB^As 



(a^-b^^-^)') 

Hence matching of the second fundamental form demands 



As fdB\ dAs {B + If As) 



(60) 



Ktt'^ = and = (iiii'e-''')^- 

In particular, if we take as the Schwarzschild metric, namely. 



(61) 



V+ : ds+^ = -A(R)dT'^ + -Tj-rdR^ + R^d^'^ 
^ ' A{R) 



with A{R) = 1 - 



then it can be seen easily by a straight forward calculation that 



(62) 



Ktt'' = and Kee^ = ( AR^] 
So the continuity of Kgg over S i.e, K00~^ — K00~ |e= simplifies to 



(63) 



2M = Rj: 



1 _ e-^^R'^ + e-2Ti?2 



(64) 



which is the definition of the mass function given in equation (11). Hence for a smooth matching of the 
interior (spherical star) and exterior space-times (Schwarzschild), the interior mass function at the surface 
must be equal to the Schwarzschild mass of the exterior vacuum. 



To study the formation of the horizon in this model we have from equation (35) 



R = m{r) 
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i.e, only one apparent horizon exists in this case. 



In particular, if the exterior geometry is Schwarzschild, then apparent horizon is the Schwarzschild horizon. 
Further, as before, the time difference between the formation of apparent horizon at comoving coordinate r and 
the time of formation of singularity of a shell at comoving coordinate r is given by (r << rs) 



(65) 

(l - 2m ^2^2^ + 5^2?; - ^b{r)v3^ 

Thus locally, if tah{r) < ts{r) then final state of collapse will be a black hole while if tah{i~) > tsi"!") then 
final singularity will be naked. The above time difference depends on the tangential pressure and total matter 
density. Thus dark energy in the form of anisotropic fluid influences the collapsing process and also the end 
state of collapse. 



tah{r) - ts{r) = - [ 
Jo 
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